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Abstract 

We study the effects of strong magnetic fields on the neutron star structure. If 
the interior field of a star is on the same order of the surface field currently observed, 
the influences of the magnetic field on the star mass and radius are negligible. If 
one assumes that the internal magnetic field can be as large as that estimated from 
the scalar virial theorem, considerable effects can be induced. The maximum mass 
of stars is arisen substantially while the central density is largely suppressed. For 
two equal-mass stars the radius of the magnetic star can be larger by about 10% ~ 
20% than the nonmagnetic star. 
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I. INTRODUCTION 



It is well known that the structure of neutron stars is mainly determined by the 
nuclear equation of state (EOS) built on the strong interactions. Recent observations 
have indicated that large magnetic fields are presented at the surface of neutron stars 
(Michel 1991; Rothschild, Kulkarni & Lingenfelter 1994; Kouveliotou et al. 1998; Woods 
et al. 1999). The dipole fields inferred from the spin-down rates of neutron star rotations 
can be up to 10 15 G (for a recent review, see Reisenegger (2001)). The strength of magnetic 
fields in the inferiority of stars remains unknown. According to the scalar virial theorem 
(Shapiro & Teukolsky 1983; Lai & Shapiro 1991) the maximum interior field strength 
could reach ~ 10 18 G for a star with R pa 10 km and M pa 1.4 M . Even larger fields may 
be expected in the core of neutron stars. Such high fields will certainly play a role when 
one evaluates the EOS of neutron-star matter. Consequently, it may causes considerable 
effects on the structure of neutron stars. 

The same problem has been addressed on the magnetic white dwarfs where a surface 
magnetic field on the order of 10 5 ~ 10 9 G and an interior field of 10 9 ~ 10 13 G are 
estimated. Both the earlier work of Ostriker & Hartwick (1968) and recent calculations 
of Suh & Mathews (2000) predicted an increase of white dwarf radii in the presence 
of internal magnetic fields. In the case of white-dwarf binary system of LB 11 146 (PG 
0945+245), if one assumes an equal mass of M — 0.9 M for the magnetic one and 
nonmagnetic one, an interior field of B pa 0.5B^ (here B e c = 4.414 x 10 13 G is the electron 
critical field) would cause the radius of the magnetic white dwarf be larger by about 10% 
than that of the normal star (Suh & Mathews 2000). 

Theoretical investigation of ideal noninteracting neutron-proton-electron (n-p-e) gas 
and interacting pure neutron matter under large magnetic fields has been carried out 
by Suh & Mathews (2001) and Vshivtsev & Serebryakova (1994). Recently, Briickner- 
Hartree-Fock calculations of spin polarized asymmetric nuclear matter have been per- 
formed (Vidana & Bombaci 2002). Based on the meson field theory several authors have 
incorporated strong magnetic fields into the equation of state of a dense n-p-e system 
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under the beta equilibrium and the charge neutrality conditions (Chakrabarty, Bandy- 
opadhyay & Pal 1997, Paper I; Broderick, Prakash & Lattimer 2000, Paper II). Evident 
changes on the EOS have been found. It was demonstrated that the nucleon anomalous 
magnetic moments (AMM) play a significant role (Paper II) which may overwhelm the 
softening of the EOS caused by Landau quantization (Paper I) to the stiffening. In the 
mean time, a dramatic increase of the proton fraction with the increase of the magnetic 
field was exhibited. However, our previous calculations (Mao et al. 2002) showed that 
with the AMM of nucleons and electrons taken into account the proton fraction was 
found to never exceed the field free case. Extremely strong fields would lead to a pure 
neutron matter rather than a proton-rich matter. In this work we will study the effects 
of large magnetic fields on the neutron star structure. We will examine the EOS of a 
dense n-p-e system with the AMM of both nucleons and electrons taken into account. 
One may argue that the electron self-energy may not change substantially in magnetic 
fields when high-order terms are taken into account. However, a systematic incorporation 
of high-order contributions beyond the AMM term is not yet clear, which will be the 
topic of our forthcoming works. Here the effects of magnetic fields on different particles 
within the considered system are treated on an equal footing. The developed EOS will 
be subsequently applied to investigate the structure of neutron stars with strong inter- 
nal fields. The paper is organized as follows: a relativistic mean-field theory approach for 
dense neutron star matter is described in Section 2. In Section 3 we present the numerical 
results. A brief summary and outlook will be finally given in Section 4. 

II. THEORETICAL FRAMEWORK 

We consider a neutron-star matter consisting of neutrons, protons and electrons inter- 
acting through the exchange of a, uj and p mesons in the presence of a uniform magnetic 
field B along the z axis. The Lagrangian density can be written as (Serot & Walecka 
1986) 
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-\o^ v + \mlu^ - • R/ 4 " + imjR M • R", (1) 

and U (a) is the self-interaction part of the scalar field (Boguta & Bodmer 1977) 

U(a) = \m 2 y + l -b{ ga af + \c{g a a)\ (2) 

In the above expressions ip and ip e are the Dirac spinors of the nucleon and electron; 
a, uj^, R^ represent the scalar meson, vector meson and vector-isovector meson field, 
respectively. A^ = (0, 0, Bx, 0) refers to a constant external magnetic field. Here the field 
tensors for the omega, rho and magnetic field are given in terms of their potentials by 

= d„u u - dvUp, (3) 
R</m/ = d/jllv — d^K^, (4) 
Ffxv = d»A v - d v A^ (5) 

and o^v = | [7^, 7„]. r is the isospin operator of the nucleon and r is its third component. 
M N and m e are the free nucleon mass and electron mass and m a , m w , m p are the masses 
of the a-, uj- and p-meson. /ijv and jiB are the nuclear magneton of nucleons and Bohr 
magneton of electrons; k p = 3.5856, K n = —3.8263 and K e = a/ir are the coefficients of 
the AMM for protons, neutrons and electrons (Greiner & Reinhardt 1994), respectively. 
The third set of parameters presented by Glendenning & Moszkowski (1991) is used as the 
nucleon coupling strengths. It gives g CT = 8.7818, g w = 8.7116, g p = 8.4635, bg^ = 27.9060, 
eg* = —14.3989. This yields a binding energy B/A = —16.3 MeV, saturation density 
Po = 0.153 fm~ 3 and bulk symmetry energy a sym = 32.5 MeV. 

Let us first consider nucleons. The Dirac equation for the nucleons in a uniform 
magnetic field can be written as 



- e^4^7^ - ^Kb/iNa^F^ - M N + g a a - g^w" - ^g^^R^ 



V> = 0. (6) 



The general solutions of the above equation in the free space, i.e., in the absence of the 
meson fields, are discussed in Appendix A. The corresponding solutions in the neutron- 
star matter can be directly written out by replacing the free quantities with the effective 



ones. The positive energy of the protons in the Fermi sea (^,5) and the negative energy 
of the protons in the Dirac sea (^,5) read as 



1 1/2 



+ g^o + ^gpRo,o > , 



1/2 



Vm* 2 + 2e5z/ + 5 A) 

(^s)_ = -| (v / m* 2 + 2e J Bz/ + 1 SA) 2 +p 
The positive-energy and negative-energy spectra of the neutrons are 
(E%) + = I [(^Pl+pl + m^ + SA) 2 + p 

(E n s )_ = - { f(^+pg + m*2 + SA) 2 + p 



ll/2 1 ] 



1/2 



1 



gw^o — 2^-^0,0 



(7) 
(8) 

(9) 
(10) 



Here A = — ^k^^x^B] S — ±1 for spin-up and spin-down particles, v is the quantum 
number of Landau levels for charged particles (Landau & Lifshitz 1977). The positive 
energy of the anti-particle is just the negative of the negative energy of the particle, i.e., 
(j% tS ) = - {K,s)_, (#s) = ~ ( E s)- ( Mao > Stocker & Greiner 1999). 

Since neutron stars are cold dense matter, we shall perform numerical calculations at 
zero temperature. The general chemical equilibrium is realized for valence particles. The 
chemical potentials of protons and neutrons are defined as 



V P = e p f + g^uo + -gpRo,o, 

/in = + g w W - -jZpR$fl- 

They are related to the respective Fermi momenta via following equations: 

fe,s) 2 = (*?) 2 - (Vm* 2 + leBv + SA) 2 , 



(11) 
(12) 

(13) 
(14) 



In the above equations, the effective nucleon mass m* = Mn — g a a. a, u Q and i?o,o are the 
mean values of the scalar field, the time-like component of the vector field and the time- 
like isospin 3-component of the vector-isovector field in neutron-star matter, respectively. 
They are obtained by solving the non-linear equations of the meson fields 



m%a + bg%a l + cgX = g a p S) 



(15) 



m o R 0fi 



gup, 
1 

- ogpPO.O- 



(16) 
(17) 



Here ps, p and po,o are the scalar density, the time-like component of the vector density 
and the time-like isospin 3-component of the vector-isovector density contributed from 
the valence nucleons, that is, from the Fermi sea. In principle, there exist additional 
contributions stemming from the Dirac sea. Here they are neglected according to the no- 
sea approximation since the renormalization of the system under the external magnetic 
field is a problem to be solved. Thus, the contributed densities are p s = Ps+Psi P = Po+Po 
and p 0i0 = Po - Po , with 
eBm* 



Ps 



Ps 



2tt 2 

m* 
4^ 



x . . . Vm* 2 + 2eBv + SA , 
> > = In 

Y , Vm* 2 + 2eBv 

Y.V n Ms -(m* + S A) 2 In 



K f,u,S + e j 



Vm* 2 + 2eBv + SA 



K f,S 



z f 



m* + SA 



Po = 



eB 
1 

2^2 



(18) 
(19) 
(20) 



s L° 



SA ( , „ ^ AN ,„ / „\2 , . m* + SA 7r.\. 

+ SA)kl s +(e]) (aresm— -))(21) 



The summation of v runs up to the largest integer for which (k P f u s) 1S positive. For 
spin- up protons v starts from 1 while for spin-down protons 0. It should be pointed 
out that here the so-called spin up and spin down are just relative notes since the wave 
functions are no more eigenfunctions of 3-component spin operator (see Appendix A), 
mainly attributed to the coupling of the spin to the magnetic field. The contributions of 
the protons and neutrons to the energy density read as 

, 2 



e = —TT 

4vr 2 ^ 



k lu,s eP f + (Vm* 2 + 2eBv + SA.) In 



arcsm 



Vm* 2 + 2eBv + SA 
m* + SA\ 7T 
J ~ 2 



(22) 



+ 



\SA -\(m* + SA) 



m* + SA)k n fS e n f + (m* + SA) 3 In 



K f,S 



(23) 



m" + SA 

In the /^-equilibrium system, the electron is assumed to move freely in the strong 
magnetic fields. The wave functions of electrons are the same as that of protons in the 
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free space except that the corresponding quantities are replaced by the electron ones. The 
energy spectrum of the electrons can be expressed as 

1/2 



El 



v,S 



ml + leBv + SAj + p 2 z 



(24) 



here A = —^K e ji B B and (^v,s) = — {pl,s) ■ The chemical potential of electrons 
li e — e e f. Its relation to the electron Fermi momentum is 



The electron density is defined as 



m 2 e + 2eBu + SA^j 



(25) 



eB 



(26) 



S v 



The contribution of the electrons to the energy density reads as 



e = ^lyy 



S v 



k)^ s e) + (f] + 2eBu + Sa) * In 



k f,u,S + e / 



ml + 2e5i/ + SA 



(27) 



Finally, we obtain the energy density contributed from the neutron-star matter 



1. 
S 1 



1 1 

+^ m l^l + ~ m l R lfl + e p + e n + e e . 



2 w u 2 

The total energy density of the system is given by 

B 2 



(28) 



8tt' 



(29) 



where the last term is the contribution from the external electromagnetic field (Landau, 
Lifshitz & Pitaevskii 1984). In the charge neutral beta-equilibrated matter, the pressure 
of the system can be expressed as (Paper II) 

B 2 



8tt' 



(30) 



The inclusion of the B 2 /8n term in Eqs. (29) and (30) is equivalent to introducing the 
electromagnetic-field tensor in the source term of the gravitational equation of general 



relativity (Bocquet et al. 1995; Bonazzola & Gourgoulhon 1996; Cardall, Prakash & 
Lattimer 2001). Numerical calculations are performed under the constraints of the charge 
neutrality pg = Po an d the /^-equilibrium p n — p p + p, e (Glendenning 1997). These two 
constraint equations together with three meson equations are solved self-consistently in 
an iteration procedure. 

III. NUMERICAL RESULTS 

Since there is no information directly available for the interior magnetic field of a star, 
we assume that the field varies from the surface to the center and adopt the following 
parametrization (Bandyopadhyay, Chakrabarty & Pal 1997) 

B(p/p ) = B surf + B cent [1 - exp (-p( p / Po y)] , (31) 

where the parameters are chosen to be (3 = 0.01 and 7 = 3. We further set B sur f = aB^, 
B cen t — ax 10 4 -B^ and take a (which should not be confused with the fine structure 
constant.) as a free parameter to check the effects of different fields. In this configuration 
the magnetic field decreases from the center to the surface of a star. We further take 
the surface magnetic field in the range of 10 13 ~ 10 15 G which is in agreement with 
the values inferred from observations. Therefore, the applied magnetic field does not 
influence the stationary configuration of a spherical neutron star considered here (note 
that the stability of an axially symmetric neutron star in large average magnetic fields has 
been discussed recently (Broderick, Prakash & Lattimer 2002; A. P. Martinez, H.P. Rojas 
& H.J.M. Cuesta 2003)). The EOS of neutron-star matter under strong magnetic fields 
is depicted in Fig. 1. It can be found that the equation of state becomes stiffer than the 
field free case at the main region of density. Stronger field leads to a stiffer EOS. For the 
very low density part the situation is just inverse. A softer EOS is obtained when the 
magnetic field is incorporated which can be more clearly seen from the upper panel of the 
figure. 

In Fig. 2 we separate the cases with and without the inclusion of the AMM effects. 
Here we have dropped the contributions from the electromagnetic field to the pressure 
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and energy density and presented the EOS for the matter part. A uniform magnetic field 
from the surface to the center of neutron stars is assumed. Dashed line denotes the case 
of B = and neglecting the AMM effects. One can see that the EOS turns out to be 
softer than the field free case at small energy densities and approaches to it in the main 
domain. The inclusion of the AMM effects here practically leads to undistinguishable 
results since at such low field the anomalous magnetic moments do not play any role. 
The whole situation changes substantially when a strong magnetic field is considered, as 
depicted in the figure by the dash-dotted line and dash-dot-dotted line for B = 5 x 10 5 -B^. 
In accordance with the findings in Paper I & Paper II the EOS becomes softer compared 
to the case of B = when the AMM is neglected and stiffer when it is included. Therefore, 
the softening of the EOS at small energy densities exhibited in Fig. 1 is mainly due to the 
fact that at the surface region we have taken a relatively small magnetic field where the 
effects of the AMM are negligible. Under the employed parametrization the field increases 
with the increase of the density. The AMM terms gradually get into the game. This leads 
to a stiffening of the EOS as shown in the region of large energy densities. 

The structure of neutron stars can be obtained by applying the developed EOS to 
solve the Tolman-Oppenheimer-Volkoff (TOV) equation for a relativistic, spherical and 
static star (Shapiro & Teukolsky 1983; Glendenning 1997). We obtain stable solutions 
for the TOV equation based on the EOS with the strong magnetic fields applied. Fig. 3 
displays the gravitational masses of neutron stars as a function of central densities, i.e., 
the sequences of stars obtained under different circumstances of interior magnetic fields 
as indicated by different curves in the figure. Rather evident effects are induced by strong 
magnetic fields. The maximum mass of stars is increased drastically in the presence of 
strong fields. For the case of a = 50 the M max is enhanced by about 40% compared to the 
field free case, which can nearly balance the effects of decreasing the maximum mass if the 
hyperon degrees of freedom are taken into account (Glendenning & Moszkowski 1991). In 
the mean time the central density is largely suppressed with the increase of fields. This 
can be seen from Fig. 1 since the equation of state becomes stiff when a strong magnetic 
field is buried in a star. Note that an enhancement of the maximum mass around 13% - 
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30% has been reported (Bocquet et al. 1995; Bonazzola & Gourgoulhon 1996; Cardall, 
Prakash & Lattimer 2001) where the authors accomplished calculations for the structure 
of axisymmetric relativistic stars. The effects of strong magnetic fields were included in the 
stress-energy tensor of the gravitational equation while neglected in the input equations 
of state. 

We show in Fig. 4 the radii of neutron stars as a function of interior magnetic field 
strength. Three cases of fixed star masses with M/M Q = 1.4, 1.6, and 1.8 are investigated. 
If the field is weak, as in the normal case, the heavier stars have smaller radii due to the 
effects of the gravitational force. However, the radius increases with the increase of the 
magnetic field strength to a large extent. When strong enough fields are presented, the 
heavier stars can even have larger radii since neutron stars with different internal magnetic 
fields may belong to different star sequences. Generally speaking, the radius of a magnetic 
star can be enhanced by about 10% ~ 20% depending on the star mass compared to the 
nonmagnetic star of equal mass. 

We have also investigated the neutron star structure under a nonvaried field strength 
from the surface to the center of neutron stars. For B = 505^: the relationship between 
the star mass and the central density is undistinguishable from that of the field free case. 
The enhancement of the star radii is found to be less than 2%. It seems to be difficult 
to observe any effects of the magnetic field by measuring the star mass and radius if 
the magnetitude of the interior field is on the same order of the surface field, though the 
influence of the surface field itself on the star properties should be pursued more closely in 
a model for matter below neutron drip (Baym, Pethick & Sutherland 1971; Lai & Shapiro 
1991). Finally, in Fig. 5 we depict the mass-radius relations of different neutron-star 
sequences with the magnetic fields considered. In order to describe the surface region a 
model for non-uniform matter at low densities (Shen 2002) should be utilized, which will 
be taken into account in our future investigations. 

IV. SUMMARY AND OUTLOOK 

Within a relativistic field theory approach we have studied the effects of strong mag- 
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netic fields on the equation of state of beta-equilibrium and charge neutrality matter 
relevant to neutron stars. The anomalous magnetic moments of both nucleons and elec- 
trons are incorporated in the model covariantly. We present the analytical expressions 
of the Dirac spinors under a uniform magnetic field. Numerical results show that if the 
magnetitude of the magnetic field is on the order of the surface field of neutron stars, the 
effects of the AMM are negligible. The EOS deviates from the field free case only at very 
low density region. If a much larger field is considered, the AMM plays a significant role 
so that the EOS becomes quite stiffer compared to the nonmagnetic case. Consequently, 
if the interior magnetic field of a neutron star is in the same level as the surface field 
inferred from pulsars, it may not cause evident influences on the star mass and radius 
as well as the maximum mass of star sequence. If ultra-strong fields do exist in neutron 
stars, considerable effects can be observed. We assume the magnetic field varies from the 
surface to the center of a neutron star. For the surface field we take the values inferred 
from pulsars. Strong fields up to 10 18 ~ 10 19 G have been considered for the center of the 
neutron star. With the employed parametrization of the field changing with the density, 
the maximum mass could be enhanced by about 40%. The central density of the star is in 
turn suppressed dramatically. For two equal-mass magnetic and nonmagnetic stars, the 
radius of magnetic one can be larger by about 10% ~ 20% than the nonmagnetic one. 

It would be interesting to check whether the enhancement of the maximum mass 
induced by strong internal fields can be balanced by the hyperon degrees of freedom 
which was known to decrease the maximum mass. Since the effects of star rotations 
increase the star radii, it is necessary to study the effects of rotations and magnetic fields 
simultaneously. Theoretically one should solve the coupled Einstein-Maxwell equations 
for axisymmetric configuration with the effects of magnetic fields taken into account both 
in the source term of the gravitational equation of general relativity and the nuclear 
equation of state. If in future astrophysical observations one can catch two pulsars with 
similar rotation periods and masses but quite different radii, a possible explanation is that 
the star having larger radius may contain a strong interior magnetic field. 

Now let us discuss several issues involved for the description of an electron in intense 
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magnetic fields. It has long been recognized that the ground-state energy of an elec- 
tron may be shifted after taking into account the AMM term. In previous calculations 
(O'Connell 1968; Chiu & Canuto 1968) one thought that the Landau levels always start 
from v — 0. Thus, the ground-state energy of an electron in the presence of the magnetic 
field can be written as 

(X 

E = m e - —fj, B B. (32) 

The vacuum becomes unstable with respect to the electron-positron pair creations when 
B > 7.6 x 10 16 G, which may have dramatic astrophysical consequences. As pointed out 
in the appendix, the Landau levels start from v — 1 for spin-up particles since the wave 
functions vanish at v — 0. It causes the ground-state energy to the form of 

E = J ml + 2eB - —fi B B. (33) 
v Zix 

The critical field for pair creations now turns out to be 2.6 x 10 20 G. However, Eq. (33) is 
valid at the level of the anomalous magnetic moment term (Jancovici 1969). Higher-order 
terms, e.g. radiative correction due to vacuum polarization effect (Schwinger 1951; 1973), 
may become effective at such large field. The incorporation of these terms in the present 
model is a problem to be addressed in future studies. 

APPENDIX A 

In this appendix we derive the Dirac spinors and energy spectra of nucleons in the 
presence of a constant uniform magnetic field along the z axis. The Dirac equation for a 
free nucleon which has an anomalous magnetic moment \i n in an external magnetic field 
can be written as 



i^d* - e^p7^ M " \wno^ - M N 



i> = 0. (Al) 



If one drops the term concerning the AMM, the solutions for neutrons are just the conven- 
tional Dirac spinors (Greiner 1990; Weinberg 1995). The Dirac theory for free electrons in 
a homogeneous magnetic field was first investigated by Rabi (1928). The wave functions 
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for charged particles without the inclusion of the AMM have been studied by several au- 
thors (Kobayashi & Sakamoto 1983; Das & Hott 1996). Johnson and Lippmann (1950) 
considered the inclusion of the AMM in the Dirac equation within a noncovariant descrip- 
tion. A covariant energy spectrum was discussed by Vshivtsev and Serebryakova (1994). 
Recently, Broderick, Prakash and Lattimer (Paper II) derived the spinors and energy 
spectra for baryons in the Dirac representation. However, their formulae are quite com- 
plicated which leads to a lack of analytical expressions for the wave functions. Here we 
resolve the problem in the chiral representation. The employed 7-matrices then become 
(Itzykson & Zuber 1980) 




ct 



75 




(A2) 
(A3) 



where cr is the Pauli matrix. Equation (Al) can be rewritten as 



—lex ■ V - e 1 + T ° a 2 Bx + (3M N + ^-^K h n N a lxv F llv 



(A4) 



In the following we define A = — ^k^hnE and consider the cases of protons and neutrons 
individually. 

Protons 



Let us specify the wave functions of protons as 



i>{X) = e- %Et+ip y y+ip > z <t> p (p y ,p z ,x) 



(A5) 



In the static system we obtain the eigenequation 

I Pz-E -M N + A 



£- -Pz-E 

-Mat + A -Pz-E 

-M N - A -f_ 



\ 



-M N - A 
Pz-E ) 



( 



\ 



0f 
I 4 4) J 



(A6) 
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Here we have defined £± = — id x =F i(p y — eBx). Provided eB > 0, one can introduce the 
eigenfunction of £+, 



wo - (£f -p 

(1/ = 0,1,2,...) 



2 V eB/ 



eB 



(A7) 



where H v (x) is the Hermite polynomial defined by 



(A8) 



Iv, Py {x) is normalized as 



y dxI llj p y (x)I fJl> p y (x) — 5 Ufl , 

oo 

H Iv, Py (x)Iu, Py (x') = 5(x - x') 

v=0 



(A9) 
(A10) 



It satisfies the following relations: 



£-4, P3 ,(z) = -W2eBuI v ^ Py {x), (I- 1:Py (x) = 0) 
i+h.pyix) = i^2eB(v + l)I u+l , Py {x). 



(AH) 
(A12) 



The eigenvalues and eigenf unctions of Eq. (A6) can be deduced in a standard way by 
performing matrix calculations. The energy spectra of protons are 



(El,s) + = E» S , 



(El,s)_ = ~E v s, 



(A13) 



with 



.1/2 



{yjM 2 N + 2eBv + SA^j +p 2 z (A14) 

here S — ±1 for the spin- up and spin-down particles. The respective eigenfunctions are 
as follows: 

-1/2 

1p!(X) = e -^ +1 t+ip y y+i Pz z ^ Kj_+Pz 1+ 'h-Bv 



\ 2E 



+i 



^JM 2 N + 2eBv + Mtv 
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X 



W2eBu 
^/Ml+2eBu+M N 



^M%+2eBv+A . . 



W2eBv ^ Mj f +2eBv+ A^j 
(E+!+P*) {^/ M 2 N +2eBu+M^j 



X 



Iu-l, Py {x) 



(A15) 



MX) 



-iE^t+ipyy+ipzZ 



1 + 



2eBv 



x 



\ 2E V _ X 

W2eBv j ( \ 

v /M^+2eBu-M N 1 ^y [X) 

^Ml+2eBv-A , , 



sjM 2 N + 2eBv - M N 
\ 



-1/2 



M%+2eBu-A 



(A16) 



■03 (X) = e - i (~ E +i) t+i Pyy +i P zZ 



Pz 



2E» +1 



1-1/2 



2eBv 



M% + 2eBv + M A 



N 
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X 



/\ 2(Bi<{ \/M'i -2. LSi'+A ') \ 



{E» +1 + Pz ) (y 'M 2 N +2eBu+M N ^j 



W2eBv 



y /Mf f +2eBu+M N 



y/M%+2eBu+A 



(A17) 



MX) 



-{(-E'^^t+ipy-y+ipzZ 



\ 2E -i 



1 + 



2eBu 



1 iV2TBl)(^/Ml+2eBu-A\ 

\ L J [ X) 

(^/M%+2eBu-M N ^j "'^ 



X 



\jMfj + 2 eJ Bz/ - M N 
\ 



-1/2 



iV2eBv 



^/M 2 N +2eBu-M N 



Iv,Py 



V 



^/M%+2eBu-A 



Pi/ 



(A18) 



/ 



One can easily check that = 1,4) forms a complete orthogonal set. Note that the 
Landau levels start at v — for spin-down particles and v — 1 for spin-up particles since 
■ipi and ^3 vanish at v — 0. As pointed out in Sect. II, here the so-called spin up and spin 
down are just relative notes for convenience of description. 

Neutrons 

The wave functions of neutrons can be specified as 



(A19) 
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Inserting it into Eq. (A4) we have the following eigenequation: 

( p z -E Px - ip y -M N + A W rhW \ 

Px + ip y -Pz-E -M N - A 

-M N + A -p z -E -(p x -ip y ) 

y -M N - A -(p x + ip y ) p z -E ) 



(A20) 



Through solving the above matrix equation we obtain the energy spectra of neutrons as 



{Eg), = E s , 



(E n s )_ = -E S , 



(A21) 



with 



Eg = 



(^pI+pI + m% + sa) 2 +p 



1/2 



The corresponding eigenfunctions read as 
MX) = e 



-iE +1 t+ip-x E +1 +p z 

V 2E +1 



1 + 



Pl+P 2 y 



( 



X 



Px-iPy 



(y/M%+p>+p} + M N y 

\ 



-1/2 



y/M 2 N +p%+pl + M N 



(S+i+Pz) ( -/Af2 +pl+^+M^J 
1 



MX) 



-iE-it+ip-x. 



I E_ 1 +Pz 

2E_! 



1 + 



pI+pI 



1-1/2 



(JW+pJ+pJ - M N ) 



(A22) 



(A23) 
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X 



Px~iPy 



^/m 2 n +pI+pI-m n 



^/m 2 n +pI+pI-a 

E-i+p z 



(Px-iPy) [^/Ml+pl+pl-I^j 
(E- 1+ p z )(^Ml+pl+pl-MN 



V 



(A24) 



/ 



MX) = 



1 + 



-i(-E+i)t+»p-x E+l+Pz 

V 2E +1 

(B+i+fc)(v^+^+?f+ M 'v 

1 



^m 2 n + pI + pI + m n ) 



-1/2 



X 



V 



y/ 'm 2 n +pI+pI + m n 



y/M 2 N +pj+pl+& 
E+i+p z 



(A25) 



MX) 



-i(-E_i)t+ip-x / -^-l +P 



1 + 



(^/M^+^+pg-M*)' 
^ (p^-ip ! /)(v /M ^+pi+pi- A ) \ 

(£- 1 +p z ) ( ^/M^ +p2 +p | _ Mjv j 



-1/2 



X 



-1 



Px-}Py_ 



^/m 2 n +pI+pI-m n 



V 



(A26) 



Again, = 1,4) is orthonormalized. 
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Figure 1: The equation of state of neutron-star matter with strong magnetic fields. The 
upper panel gives the relationship between the pressure and total energy density while 
the lower panel displays the pressure as a function of density. Different curves correspond 
to different magnetic field strength as indicated in the figure and explained in the text. 
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Figure 2: The equation of state of neutron-star matter with a uniform magnetic field. 
The contributions of the magnetic field to the pressure and energy density are neglected. 
Different curves are related to different cases of magnetic fields B and with or without 
the inclusion of the anomalous magnetic moments k as indicated in the figure. 
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Figure 3: The gravitational masses of neutron stars as a function of central densities, 
i.e., the sequences of stars under different interior magnetic fields as depicted by different 
curves in the figure. 
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Figure 4: The radii of neutron stars with equal gravitational masses as a function of 
interior magnetic field strength. Different curves correspond to different groups of star 
masses. 
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Figure 5: The mass-radius relations of neutron-star sequences. The different curves are 
related to the different interior magnetic fields as depicted in the legend of Fig. 3. 
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